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Abstract 

The N —* oo limit of the edges of finite planar electron densities is discussed for 
higher Landau levels. For full filling, the particle number is correlated with the mag- 
netic flux, and hence with the boundary location, making the N — > oo limit more 
subtle at the edges than in the bulk. In the n th Landau level, the density exhibits n 
distinct steps at the edge, in both circular and rectangular samples. The boundary 
characteristics for individual Landau levels, and for successively filled Landau levels, 
are computed in an asymptotic expansion. 



1 Introduction 



Recent work ][], 0, f|, [j[ || [7j has illustrated the fundamental importance of the edge proper- 
ties of planar incompressible quantum fluids in the understanding of the quantum Hall effect 
||. The physical incompressibility of the quantum Hall samples correlates boundary and 
bulk degrees of freedom in an interesting manner. The boundary values of the bulk particle 
densities behave as 1 + 1 dimensional chiral Kac-Moody currents, the algebraic properties 
of which characterize the edge excitations P, E3, EL Furthermore, the incompressibility may 
be fruitfully understood via a quantum deformation of the algebra of area-preserving dif- 
feomorphisms |6], [7], [|. Most of this work has concentrated on the physics of the lowest 
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Landau level, but higher Landau levels are also important for a complete description of the 
quantum Hall effect. Here, I consider the exact large N asymptotics of the edge properties 
of incompressible quantum fluids in higher Landau levels. 

Consider the nonrelativistic Pauli hamiltonian operator for spin-polarized (i.e. scalar) 
fermions confined to the two-dimensional plane, in the presence of a uniform magnetic field 
(of strength B) perpendicular to the plane: 

H=^(p--A) -— - (1) 



2m \ c ) 2mc 

It is well known that, in the plane, the spectrum of this hamiltonian consists of an infinite 
discrete set of equally spaced energy levels (called "Landau levels"), each of which is infinitely 
degenerate []10| |. The infinite degeneracy of each level is directly related to the infinite area 



of the plane, and for a sample of finite area the degeneracy N is of the order of the net 
magnetic flux through the region (measured in fundamental flux units hc/e): 

iV = + /i (2) 

where /i = 0(1). In the plane, the eigenstates of the energy operator (|I|) may be written as 
ipi, where n = 0,1,2,... is the Landau level index (corresponding to the energy eigenvalue 
E n = rihuj c where u> c = eB/mc is the cyclotron frequency) and j labels the degenerate states 
within each Landau level. In the n th Landau level, the expectation value of the density 
operator is given by the function 

Pn(x)=Y,\r n (x)\ 2 (3) 

3 

For a fully filled Landau level, the occupation number equals the degeneracy, so the sum- 
mation in (|3|) is always over N values of j. At some fixed point x within the bulk of the 
sample, the thermodynamic limit, N — > oo, simply involves extending this finite sum to an 
infinite sum. But at the edge of the sample, each individual density \ip J n {x)\ 2 in the sum is 
itself evaluated at a distance scale \x\ = R which is correlated with the particle number as 
R ~ ~ yN because of (fj) and the condition of full filling. Thus the thermodynamic 
limit is much more subtle at the sample boundaries. The analysis of these limits is the 
subject of this paper. 

The precise form of the eigenstates ij)° n depends on the gauge chosen for the vector po- 
tential A. To describe a finite sample, one should choose the gauge with the geometry of 
the sample boundary in mind. It is convenient to choose the gauge so that, on the bound- 
ary, the normal component of the vector potential A vanishes. Thus, for circular samples 
one uses the "symmetric" gauge, A = —^(y,—x), while for rectangular samples (actually 
for rectangular strips with upper and lower sides identified) one uses the "Landau" gauge, 
A = B(0, x). Owing to the degeneracy of the Landau levels, a change in the gauge involves 
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a change in the nature of the sum over j in ([3]), and so the functional form of p n is gauge 
dependent. We shall see that it is considerably easier to analyze the asymptotics of p n in 
the symmetric gauge than in the Landau gauge. 



2 Circular Geometry 

In the symmetric gauge it is useful to work with complex coordinates z = x + iy, and the 
degenerate states are labelled by an (integer) angular momentum index j: 

In (HD, the Landau level index n runs over all nonnegative integers n — 0, 1, 2, . . ., while the 
angular momentum index takes values j = —n, —n+1, —n + 2, . . . within the n th energy level. 
The L J n are the generalized Laguerre polynomials |Tl| and the eigenstates ip 3 n are orthonormal 



in the plane. The individual densities \ip^(x)\ 2 are radially symmetric functions of Gaussian 
form, peaked at a radius r ~ y/j + n, which increases with j. So for a finite disc [] of TV 
electrons, all in the n th Landau level, the density function is 



N-n-l 



Pn (t,N) = n\ e* £ jJL- (m)T (5) 

where the natural dimensionless variable is ^ = -|r 2 , in terms of which the boundary r = R 
is given by £ = -f -R 2 = 0. Thus, analysis of the boundary density in the thermodynamic 
limit involves the iV — > oo behaviour of p n (£, A/") for ^ = <p ~ AT. 

Consider first of all the density function for the lowest Landau level, 

2V-1 p 

Po(e,iV) = e^E^ 

3=0 

= e^e^iCO, (6) 
where ejv-i(0 is the truncated exponential function ffTTH . This lowest Landau level density 



function is plotted in Figure ([I] a) for A^ = 100. We see that po(£> A/") does indeed represent 
a finite droplet of uniform bulk density (equal to 1 in the normalization chosen in equation 
(P)), dropping rapidly to zero at the 'edge' where £ ~ N. For fixed £ within the bulk it is 
clear from @ that po(£> AT) — > 1 as A^ — > oo. Indeed, since the bulk density is asymptotically 
uniform (see below), we can evaluate this uniform value deep in the bulk (i.e. at £ = 0) and 
(f|) immediately gives po(0, N) = 1. 

a similar way, one can of course treat an annulus as the difference of two discs. 
2 The overall normalization has been chosen so that J d£p n — N, the degeneracy number. 
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To study po(£j N) in the vicinity of the boundary (i.e. evaluated at £ = = iV — p), we 
use the following integral representation | ij for the truncated exponential function in terms 
of the incomplete Gamma function T(N, £) 

ev-i(0 = -Y^e^ (7) 
where _ 

poo 

r(N,£)= / e-H^dt (8) 
One can then use Laplace's method |L2| to show that the large N behaviour near the edge is 

giving a leading order edge density equal to half the bulk density, with subleading corrections 
that go as 0(-^=) in the variable £. The £ derviative of the lowest Landau level density is 

= - (iV^T)! e (10) 

Note that the summation over j has disappeared, due to cancellations between successive 
terms, a fact that makes the iV — > oo limit of the derivative much easier to consider. For 
large N, the derivative is nonzero only in the vicinity of the boundary. In fact, for large N, 
all derivatives of po vanish in the bulk, explaining the uniform bulk density. We can estimate 
the derivative a t the boundary using the fact 

E N * 1 (g-iv) 2 / 1 \ / 1 (£ - AH 3 \ 

Therefore, 



N\ y^N V 12iV J \ 3 A^ 2 



^-^-^ l2N { 1 + ^ + -) (12) 



However, since the sample itself grows in size as iV increases (R ~ v AT), it is more meaningful 
to consider the derivative in terms of the rescaled variable ( = £/N, for which the boundary 
is at £ = (1 — jr). In terms of this rescaled coordinate, 

d :Po\^ /N = -M^- NWN)2/2 ( 1 + 0(h) (13) 



which shows that the edge derivative becomes large and negative as —\<N; and that the 
derivative is a delta function of width concentrated at the boundary, as is appropriate 
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for the derivative of a step-function density profile. In terms of the actual radial coordinate 
r = \/-§£, the slope at the boundary is asymptotically constant 

d 

dr' 




po| boundary = -y- i + (i4) 



These results for the lowest Landau level density are reassuring, but not particularly sur- 
prising. The N — > oo limit is greatly facilitated by the relation of po(£> N) to the incomplete 
gamma function. Such a direct relation is not available for the higher Landau levels due 
to the appearance of the Laguerre polynomial factors in (|5|). However, the asymptotic edge 
behaviour is significantly more interesting for these higher Landau levels. In Figure ([I]), 
the densities p n are plotted for n = 0, 1,2,3 with iV = 100. The pattern is striking - each 
density function has uniform bulk value equal to unity, and each falls off at an edge located 
at £ ~ N, but the n > 1 densities exhibit a distinct step-like pattern. This is a general 
feature of p n (£, N) - there are exactly n such steps in the vicinity of the boundary. However, 
if one looks more closely by magnifying the scale of p n considerably, these " steps" are in fact 
pairs of a local minimum and a local maximum. The locations of these local minima and 
maxima may be identified as a result of the following remarkable identity which generalizes 
the formula fllPp: 

J d-JV-n-l 

-p n (£, N) = -^f^YyLn^iO L N n - n {i) (15) 

Note that, as in (|10D, taking the derivative of p n enables the removal of the summation 
over the angular momentum index: this is due to intricate cancellations between successive 
terms in the sum, using the recurrence and differential-difference relations of the Laguerre 
polynomials. From (|15|) we see that the 'troughs' and 'peaks' of each step are located at the 
zeros of L^~ n+1 (^) and L^~ n (£). Each Laguerre factor in ( |T5D is a polynomial of degree n 
with n real zeros, and, furthermore, their zeros are interwoven (and for a given n, the k th 
zero of JJ n increases with j |]TT|), explaining the pairing of local maxima and minima which 



gives the impression of steps. 

Perhaps even more remarkable than this step-like behaviour is the fact that when the 
density functions for successive Landau levels are summed, the boundary again becomes 
"smooth" (i.e. the steps disappear). This is directly relevant for applications to the quantum 
Hall effect ||, where one is interested, for example, in the density for full filling of the first 
p Landau levels: 

P%Lit,N) = T,Pn{t,N) (16) 

These fully filled densities are plotted in Figure (§) for p = 1, 2, 3, 4 and for iV = 100. Notice 
that the form of these density functions is once again like the lowest Landau level density 
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function po plotted in Figure a). To explain this miraculous cancellation of the boundary 
steps, note first of all the following identity 3 which relates the density in the n th level to the 
density in the (n — l) th level: 

tN-n 

p n (t,N)=p^ l {Z,N) - (n-l)\e-Z^—y L»- n (OL^r(0 (17) 

The second term on the right hand side is concentrated at the boundary, since it is a sum 
of terms like (0). Using this identity, it follows that 

/&L(£, N)=P M£, AO - |> - n)(n - l)!^-"^-"(0^-T(0 (18) 

Once again, the sum on the right hand side is concentrated at the boundary, while the bulk 
value is just p (i.e., p times the bulk value of p ). For low values of p it is easy to perform 
the finite sums in (|T8|), yielding 



tN-l 

f%L(Z,N) = 2 Po (e,iV) + e-^ ( ^_ i)! (e-iV) 

cN—2 

pSLfoN) = 3 Po (e, AO + e~* ^ _ (e - (3N - 5)£ 2 + N(3N - 7)£ - iV(iV - 1) 

pSaiCCiV) = 4p (e, iV) + _ - (5iV - 8)e 4 + (10iV 2 - 30N + 26)£ 3 

-iV(10iV 2 - 42iV + 50)£ 2 + N(5N 3 - 26iV 2 + 37iV - 16)f 
-N(N -1) 2 (N -2) 2 ) (19) 

The edge contributions are all concentrated at the boundary, and despite the seemingly 
complicated expressions (|i~9|), they all have a common functional form - see Figure ([!]). 

Thus, the leading behaviour of the density of p fully filled Landau levels is identical 
to p times the density of the fully filled lowest Landau level, with subleading corrections 
~ 0(1/ \/~N) concentrated at the boundary. For full filling of the first few Landau levels, this 
subleading edge behaviour can be determined from (Q) and (|T9"|): 



3 This is a nontrivial identity which relies on the recurrence relation properties of the Laguerre polynomials. 
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The general formula for this subleading behaviour, for any p, is 
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(20) 



(21) 



where [2^-] denotes the integer part of Notice that the coefficient of j!^ N for p odd is 
equal to that for p + 1. 



The derivative of PtoLi * s localized at the boundary, since 



d_ 
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Ptotal 
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(22) 



For low values of p, this derivative is: 
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+4(iV - 1)(N - 2)(5iV - 3)£ 3 - 3N(N - l)(iV - 2)(57V - 11)£ 2 
+QN(N - l) 2 (iV - 2)(N - 3)f - iV(iV - l) 2 (iV - 2) 2 (iV - 3)) (23) 
In terms of the radial coordinate, the asymptotic boundary slopes are 
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The general formula for the leading-order boundary slope is 



(p) I boundary 

dr Ptotaxl 



r B2r( 
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(25) 



where, as before, f 2 ^-] denotes the integer part of Notice that for p odd the asymptotic 
boundary slope is equal to that for p + 1. Note also that the leading dependence of the slope 
agrees with the 0(p/ \/N) dependence of p^ ta \(N — p,, N) given in (^[Ti]), as it must. 

These results rely on such specific properties of the Laguerre polynomials, that it is not 
immediately clear whether these strange edge phenomena are some sort of side-effect of 
the symmetric gauge or of the circular nature of the samples. It is therefore instructive to 
consider also the boundary phenomena for the rectangular samples using the Landau gauge. 
Notice that the Landau and symmetric gauges are related by the gauge transformation 
^Landau _ ^symmetric _|_ with gauge function A = y xy. Normally, in the absence of 
degeneracies, different gauges lead trivially to identical densities since the nondegenerate 
wavefunctions are simply related by a phase factor e tX . However, due to the degeneracy 
of the Landau levels, the gauge invariance only implies that e tX ^ s y mmetric j s some linear 
combination of the ^ Landau with the same energy (i.e. with the same Landau level index 
n). The particular form of this linear combination is gauge dependent, and so the conversion 
of the density from one gauge to another is not such a straightforward matter. 



3 Rectangular Geometry 

In the Landau gauge, the eigenstates of the Pauli energy operator flU) are 

€ = e ikV f-V /4 ~4^H n (VB(x - 4)1 exp(-^[x-^\ 1 (20) 



7i-y 




n! 

where, as in (f|), the Landau level label n = 0, 1,2, . . ., but now the degenerate states are 
labelled by a continuous momentum index k. This continuous momentum label may be 
discretized by compactifying in the y direction. That is, if we consider the infinite strip 
— oo < x < +oo, < y < L rather than the entire plane, then with periodic or anti-periodic 
boundary conditions in the y direction, the momentum index k takes values k = (^)j- 
Here, j takes integer values for periodic boundary conditions, and half-odd- integer values 
for anti-periodic boundary conditions. The H n are the Hermite polynomials |TT |, and the 
wavefunctions ipn are orthogonal in the infinite strip. 

The individual densities \4>n\ 2 are independent of the y coordinate (just as in the sym- 
metric gauge the individual densities are independent of the angular coordinate), and have 
a Gaussian form in the x direction (analogous to the radial Gaussian form of the densities 
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in the symmetric gauge) peaked at x ~ -g. Thus, for a rectangular strip of finite extent in 
the x directionfj there are corresponding upper and lower bounds on the allowed values of 
the discrete momentum.^] For convenience of notation, we choose the sample to be a square: 
— ^ < x < < y < L, and use the rescaled x-coordinate, £ = ^-x = jx, in terms of 
which the boundaries lie at £ = ±|. 

Then the density for the fully filled n th Landau level is 




Pn&N) = J-— Y e - 2 ^)/^ H n L -(£-j) (27) 




where the flux and the particle number N are related as in (Q) for full filling. 

The density ( f2"T| ) is plotted in Figure (|J) for the first four Landau levels for N = 40 and 
~ iV.[] From these plots we observe the same behaviour as in the symmetric gauge. In the 
lowest Landau level, there is a uniform bulk density equal to 1, which drops rapidly to zero 
in the vicinity of the boundaries at£ = ±|~±y. In the higher Landau levels there is once 
again distinct step-like behaviour at the edge - specifically, in the n th Landau level there 
are n distinct steps in the density in the vicinity of the boundaries. Moreover, when the 



densities for successive Landau levels are summed, as in (16), these boundary steps cancel 
out, leaving a smooth boundary density of the same form as the lowest Landau level case: 
see Figure (|5]). 

It is much more difficult to prove rigorous results for the asymptotics of the Landau gauge 



p n (£,iV) because there are no clear analogues of the Laguerre polynomial results ( p!5| , p^ ) 
for the Hermite polynomials. However, the Christoffel-Darboux identity for the Hermite 
polynomials JTT[ implies the following finite sum formula: 



= 2^! ((^W) 2 " BmWH^xj) (28) 

Applying this formula to the fully filled density PtoLi = 2=o Pn with p n as in (|27|) we see 
immediately that 

N / 2 

(p) 

Ktotal 



(29) 



4 Note that this means that the topological equivalent is an annulus rather than a disc, but an annular 
density may be simply treated as the "difference" of two discs. 

5 This is Landau's original argument for estimating the degeneracy in terms of the magnetic flux through 
the finite sample JlO[ ] . 

6 Taking <f> = N — [i with [i = 0(1) does not affect these plots in any noticeable manner. Also note that I 
have taken N — 40 rather than N = 100 simply due to a limited graphics-plotting capacity. 
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The summation term on the right hand side of (|29| ) is concentrated on the boundary - the fact 
that it vanishes in the bulk can be seen most directly by noting that the orthogonality of the 
polynomials H p _i and H p+ % ensures that the leading contribution to the Euler-MacLaurin 
formula |12[ vanishes. Therefore, the bulk value of the fully filled density is asymptotically 
equal to p times the bulk density of the p th Landau level. This bulk value (for any given 
Landau level) may be evaluated deep in the bulk at £ = by applying the Euler-MacLaurin 
formula to p n (0, N) 



'2 I l r+ N / 2 



1 + 



N , 



(30) 



where we have used the fact that ^==y |e x2 / 2 H n (x)\ < k for all x and n, where k = 1.086435... 
is Charlier's number [lT|. By a similar computation, now evaluating p n {£,N) at the edge 
where £ ~ ±y, we find 



N 




+ 




(31) 



4 Conclusion 



In this paper it has been shown that, in the iV — > oo limit, the expectation value of the planar 
electron density for the n th Landau level exhibits n distinct steps near the boundary. This 
has been shown in both the symmetric gauge (appropriate for circular samples) and in the 
Landau gauge (appropriate for rectangular samples), in order to dispel possible doubts that 
this is perhaps an artifact of the particular (gauge dependent) basis used for the degenerate 
subspaces. At first sight, such non-smooth boundary structure seems to be problematic for 
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the usual |2|, |3], [|, [!]] identification of the boundary charge with a chiral Kac-Moody current. 
However, when successive Landau levels are filled these boundary steps cancel out, leaving a 
smooth boundary whose boundary characteristics have been computed exactly (see equations 
(|2l|j25|) ). It is not a priori obvious, in this independent particle picture, that the densities 
for independent Landau levels must conspire to produce a smooth edge behaviour. 

Having identified the boundary behaviour of the higher Landau level electron densities, 
it would now be interesting to consider the boundary currents within a given Landau level, 
as discussed, for example, in || for the lowest Landau level. Presumably, for successively 
filled Landau levels one obtains different representations of the Wi +OQ algebra. It would be 
very interesting to learn whether there is any connection between the boundary steps of the 
n th Landau level density and the different "boundaries" needed for the Wi +OQ representation 
theory and for the fractional QHE hierarchy [|, ^|. It should also be possible to apply 
some of the results obtained here, for the asymptotic behaviour of sums of Laguerre and 
Hermite polynomials, to the analysis of the pair correlation function |T3| and to the Coulomb 
interaction energy in higher Landau levels. 
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Figure 1: The symmetric gauge density functions (|5|) plotted for N = 100 electrons for 
the first four Landau levels, n = 0, 1,2,3. Notice the uniform bulk density equal to 1 and 
the boundary fall-off at £ ~ N. For the n th Landau level there are n distinct steps at the 
boundary. The exact locations of these steps are given (see equation ([15])) by the zeros of 
Ln~ n ~ 1 and L%~ n . 
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Figure 2: The symmetric gauge density functions Ptotai(£> N) for p fully filled Landau levels, 
plotted for N = 100 electrons. Notice the uniform bulk density equal to p, the boundary 
density equal to |, and the absence of the boundary steps seen in the individual higher 
Landau level densities in Figure ([[]). 
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Figure 3: The N = 100 edge contribution by which the symmetric gauge density for p fully 
filled Landau levels differs from p times the density for the lowest Landau level - see ( JT8| ; [T9| ) . 
Notice that the edge contribution is centered at £ ~ 100 and that for different p it has the 
same functional form. 
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Figure 4: The Landau gauge density functions ( p7|) plotted for N = 40 electrons for the 
first four Landau levels, n — 0, 1,2,3. Notice the uniform bulk density equal to 1 and the 
boundary fall-off at £ ~ ±y. For the n th Landau level there are n distinct steps at each 
boundary. 
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Figure 5: The Landau gauge density functions Ptotai(£> -^0 f° r P f un Y filled Landau levels, 
plotted for N = 40 electrons. Notice the uniform bulk density equal to p, the boundary 
density equal to |, and the absence of the boundary steps seen in the individual higher 
Landau level densities in Figure (|). 
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